Inflationary spectra and partially decohered distributions 



in 
o 
o 

(N 
-i— > 

o 
O 



(N 

> 



CO 

c3 



David Campo* 
Department of Applied Mathematics, University of Waterloo, 
University Avenue, Waterloo, Ontario, N2L 3G1 Canada 

Renaud Parentani^ 

Laboratoire de Physique Theorique, CNRS UMR 8627, 
Bdtiment 210, Universite Paris XI, 91405 Orsay Cedex, France 



Abstract 



It is generally expected that decoherence processes will erase the quantum properties 
of the inflationary primordial spectra. However, given the weakness of gravitational inter- 
actions, one might end up with a distribution which is only partially decohered. Below a 
certain critical change, we show that the inflationary distribution retains quantum prop- 
erties. We identify four of these: a squeezed spread in some direction of phase space, 
non-vanishing off-diagonal matrix elements, and two properties used in quantum optics 
called non-P-representability and non-separability. The last two are necessary conditions 
to violate Bell's inequalities. The critical value above which all these properties are lost 
is associated to the 'grain' of coherent states. The corresponding value of the entropy 
is equal to half the maximal (thermal) value. Moreover it coincides with the entropy 
of the effective distribution obtained by neglecting the decaying modes. By considering 
backreaction effects, we also provide an upper bound for this entropy at the onset of the 
^ . radiation dominated era. 



INTRODUCTION 



X 

Inflation tells us that the primordial density fluctuations arise from the amplification 
of vacuum fluctuations [1—5]. As a result of this amplification, the initial vacuum state 
becomes a product of highly squeezed two-mode states [6]. In spite of their intricate 
character, expectation values are rather simple. For instance, when evaluated on the 
Last Scattering Surface, besides the small deviations from a scale invariant spectrum [7], 
the non-trivial information delivered by the two-point function concerns the temporal 
coherence of the modes [8]. This last property follows from the fact that its Fourier 
transform, i.e. the power spectrum, displays oscillations and zeros. Since we are dealing 
with an ensemble, the presence of these zeros tells us that all realizations of the ensemble 
with a given wavenumber have the same temporal phase. In classical terms, this coherence 
can be easily enforced by neglecting the decaying mode. The residual random properties 
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of the distribution then only concern the amplitude of the growing modes. In the high 
occupation number limit, these amplitudes can be treated as stochastic variables. 

In brief, when considering the physics which took place near the recombination, it is 
convenient and sufficient to neglect the decaying mode. However this simplified description 
has several drawbacks. Indeed, the settings are too restrictive to describe distributions 
which result from decoherence processes (or more generally any non-linear process) taking 
place in the early universe. In particular, generic modifications of two-point function from 
the standard result cannot be parameterized in these simplified terms. This is true even 
for modifications which preserve the isotropy and the Gaussianity of the distribution. 

The aim of this paper is to analyze quantum distributions characterized by an arbitrary 
level of coherence. We identify the parameter which governs this level and determine how 
the properties of the distribution depend on it. At fixed occupation number, two-mode 
squeezed states give the most coherent distribution, and thermal states the least one. 
Starting with squeezed states in the high occupation number limit, when slightly modifying 
this parameter, the properties are either extremely robust, and therefore accessible in the 
classical regime, or extremely fragile. One finds that the fragile properties are quantum 
mechanical in character. They concern the squeezed spread in some direction in phase 
space (also called the sub-fluctuant direction), the quantum coherence of macroscopically 
different states in the perpendicular direction, as well as two refined properties called the 
non-separability [9] and the non-P-representability [10] of the distribution. We show that 
both the last properties are necessary to have expectation values which can violate Bell's 
inequalities. Below a critical value of the decoherence parameter, the distributions retain 
quantum properties, and can therefore violate these inequalities. 

Two-mode coherent states provide an appropriate basis to perform this analysis. The 
reasons are the following. First, as well known, a coherent state provides the quantum 
counterpart of a particular classical realization of an ensemble. This correspondence is 
well defined for the highly excited (amplified) modes we are dealing with, because the 
spread in position and momentum is negligible with respect to the "displacement" in phase 
space which encodes the occupation number. Remember that the observed temperature 
anisotropies of relative amplitude 1CT 5 fix the mean occupation number to be of the order 
of 10 100 . 

Second, homogeneity and isotropy restrict the non-vanishing matrix elements of the 
distribution to two-mode elements with opposite wave vectors k. Hence, superpositions of 
two-mode coherent states can characterize Gaussian, homogeneous, and isotropic distri- 
bution with arbitrary level of coherence. Two important relationships should be noticed. 
First, the level of coherence is fixed by the expectation value of two destruction operators 
with opposite wave vectors. (This expectation value determines the number of quanta 
in entangled pairs.) Second, the residual coherence of the distribution is related to the 
residual value of the decaying mode. Therefore the 'quality' of temporal coherence of the 
modes at horizon re-entry can be (in principle) traced back to the quantum distribution. 

Third, two-mode coherent states allow to make contact with the general remark of 
[11, 12] according to which squeezed states rapidly decay into a statistical mixture of 
coherent states when the oscillators are weekly coupled to an environment. When no longer 
neglected, small non-linearities [13] amongst modes in inflationary models will also induce 
some decoherence. (Let us make clear that the apparent violation of unitarity follows, 
as usual, from the neglect of the correlations among modes when evaluating the entropy, 
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see footnote p. 188 in [14]. If one keeps track of these non-linearities, the non-linear 
evolution would give no entropy increase.) When applying the considerations of [11, 12] 
to two-mode squeezed states, the distribution becomes diagonal in the basis of two-mode 
coherent states, but not of one-mode coherent states. As we shall see, this guarantees 
that the temporal coherence of modes is preserved while the quantum properties are 
progressively erased. 

The distribution which is diagonal in the two-mode coherent states basis plays a critical 
role because it separates distributions which have kept their quantum properties from 
those which have lost them. The corresponding value of the entropy is 1/2 the maximal 
(thermal) value, given the occupation number. 

To establish the critical character of this decoherence scheme we explore the space of 
(partially) decohered Gaussian distributions. When dealing with macroscopic occupation 
numbers, the properties of these distributions fall into two separate classes: they are ei- 
ther extremely robust or extremely sensitive to the level of coherence. (In more dynamical 
terms, this would translate into a insensitivity/sensitivity to a weak coupling to an envi- 
ronment.) The entropy is found to be a sensitive variable whereas the power spectrum is 
a robust one: the relative modifications of the latter are 1/n whereas the changes of the 
former are in Inn. (Throughout the paper we shall evaluate these entropies exactly by 
exploiting the fact that all Gaussian distributions can be expressed in terms of thermal 
distributions of new oscillators, see Appendix C.) Similarly, the residual squeezing and 
the quantum coherence are both sensitive whereas the temporal coherence of the modes 
is the second robust property of the distribution. In fact, we shall see that all quantum 
properties are lost for a critical value of the residual coherence of the distribution. This 
critical value is that given by the grain of coherent states. 

We then exploit the sensitivity of the entropy to establish a correspondence between 
the above distribution and the effective one obtained by neglecting the decaying mode 
and used in numerical codes such as CMBFAST [15]. Since their entropy coincide, and 
since the entropy faithfully traces the residual coherence, this agreement tells us the quan- 
tum distributions which correspond to the effective one are very 'close' to the diagonal 
distribution in the two-mode coherent state basis. 

Finally, we also provide an upper bound for the entropy by considering backreaction 
effects at the end of the inflationary period. This upper bound is given by 3/4 of the 
maximal entropy. 

Having established these properties, two questions should be confronted. The first 
concerns the efficiency of decoherence processes which occurred in nature: Were they 
powerful enough to reach the critical decoherence level associated with coherent states ? 
This question is currently under investigation. Preliminary results indicate that the critical 
value is indeed reached thereby implying that the resulting entropy is larger than the above 
bound. The second question is more important. It concerns the impact of the decoherence 
level on structure formation (i.e. non-linear evolution). This is a difficult question but 
we can already acknowledge the fact that the critical decoherence level associated with 
the loss of quantum coherence is so small that it could hardly have any impact on this 
evolution. Therefore the question of the existence of residual quantum properties seems 
to be of pure academic delight. 

We conclude this introduction by mentioning that closely related questions have been 
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already analyzed in several papers, see [16-23]. What we add 1 is a further clarification of 
the matters, the critical role played by coherent states in separating density matrices, the 
upper bound of the entropy, and the algebraic relationships between the entropy and the 
other quantum properties. 



II. REVIEW OF THE STANDARD DERIVATION OF PRIMORDIAL SPECTRA 
A. Quantum distribution of two-mode states 

In this subsection we recall how the amplification of modes of a massless field propa- 
gating in a FRW universe translates in quantum settings in the fact that the initial ground 
state evolves into a product of highly squeezed two-mode states. Before proceeding, we 
remind the reader that it has been shown that the evolution of linearized cosmological 
perturbations (metric and density perturbations) reduces to the propagation of massless 
scalar fields in a FRW spacetime [24]. In this article, we shall only consider the mass- 
less scalar test field since the transposition of the results to physical fields represents no 
difficulty. Indeed, when preserving the linearity of the evolution, the only modification 
concerns the late time dependence of the modes [25]. 

Let us work in a flat FRW universe. The line element is: 

ds 2 = a{ V ) 2 [-dr] 2 + Sijdx^x^ . (1) 

For definiteness and simplicity, we consider a cosmological evolution which starts with an 
inflationary de Sitter phase and ends with a radiation dominated period. When using the 
conformal time r\ to parameterize the evolution, the scale factor is respectively given by 

aM = - H( V -2 Vr ) ' fOT "°°<^<^. ( 2a ) 
a(v) = ^2 for V>Vr, (2b) 

where r\ r > designates the end of inflation. The transition is such that the scale factor 
and the Hubble parameter are continuous functions. This approximation based on an 
instantaneous transition is perfectly justified for modes relevant for CMB physics. Indeed, 
their wave vector k obeys krj r ~ 10~ 25 ~ e~ 60 when inflation lasts for 60 e-folds. Hence, 
the phase shift they would accumulate in a more realistic smoothed out transition would 
be completely negligible. 

Let £(77, x) be a massless test scalar field propagating in this background metric. It is 
convenient to work with the rescaled field 4> = a£ and to decompose it into Fourier modes 

/ ikx 
d3k j^^Mv)- (3) 



1 The present work is based on our former paper [23]. For completeness and to ease the reading, we have 
included some of its material and shall no longer refer to it. 
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The time dependent mode </>k obeys 

^ k +^ 2 -^^ k = 0, (4) 

where k = |k| is the norm of the conformal wave vector. 

In our background solution, k 2 — d 2 a/a is negative during the de Sitter period when 
the wavelength is larger that the Hubble radius. This leads to a large amplification of 4>k- 
In quantum terms this mode amplification translates into spontaneous pair production 
characterized by correspondingly large occupation numbers. 

To obtain the final distribution of particles, one should introduce two sets of positive 
frequency solutions of Eq. (4). The in modes are defined at asymptotic early time, and 
the out ones at late time. Both have unit positive Wronskian in conformity with the usual 
particle interpretation [26]. One gets 

* w - Wk i 1 ~ sorW) ■ for " K * ■ <5a> 

#r'(>)) = 4 f ^ i ", for X>.fr. (5b) 



In spite of the time dependence of the background, these two positive frequency modes are 
unambiguously defined (up to an arbitrary constant phase which drops in all expectation 
values and which has here been chosen so as to simplify the forthcoming expressions). In 
the de Sitter epoch, there is no ambiguity for relevant modes if inflation lasts more than 
70 e- folds, see [27] for the evaluation of the small corrections one obtains when imposing 
positive frequency at some finite early time. Moreover, there is no ambiguity for the initial 
state of relevant modes: at the onset of inflation they must be in their ground state [24, 28]. 

In the radiation dominated era there is no ambiguity either because the conformal 
frequency is constant since d 2 a = 0. In this era, the decaying and the growing modes 
correspond, respectively, to the real and the imaginary part of the out modes Eq. (5b), 
as can be seen by considering the limit krj <C 1. (When considering only the leading order 
terms (in l/kr] r ) of expectation values, several definitions of the growing modes can be 
used and lead to the same results. However in order to obtain unambiguous answers to 
the questions we shall ask (which concern next to leading order terms), we must adopt 
a precise definition. In the sequel, we define the growing mode by the imaginary part of 
Such a definition is necessary to relate, for instance, the power of the decaying mode 
to that of the so-called subfluctuant mode.) 

The in and out modes are related by a Bogoliubov transformation 

W(r l ) = a k cj>l ut (r l ) + m ut *(ri), (6) 
where the Bogoliubov coefficients are given by the Wronskians 

a fc = (^,C), fl = - W.C) • CO 
These overlaps should be evaluated at transition time r] r since modes satisfy different 
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equations in each era. One gets 

e 2ikr) r _ 

«=2«ir (8b » 

Thus, for relevant modes, krf r ~ 10~ 25 , the in modes are enormously amplified. Concomi- 
tantly, they are dominated by the sine during the radiation dominated era because the 
phase of — P%/otk ls much smaller than one: 



,,■„. , i [sin/cry .o cos/cry 

4>k(v) = -PT^ ~F=- + 0(kr) r ) —=J- 



V>Vr- (9) 



To leading order, i.e. neglecting the cosine in Eq. (9), the physical modes 4> % k n /a correspond 
to the growing modes. They possess a well-defined temporal behavior, e.g. they are 
constant until they start oscillating as they re-enter the Hubble radius when krj ~ 1. 

Lets now see how these considerations translate in second quantized settings. Each 
mode operator is decomposed twice 

Mv) = aitiiv) + alUnv) , (10) 

where j stands for both the in and out basis. The operators so defined are related by the 
transformation 

a£ = & &3 f . (ii) 

This transformation couples k to — k only. Hence, when starting from the in vacuum (the 
state annihilated by the a^ 1 operators), every out particle of momentum k will be accompa- 
nied by a partner of momentum — k. Moreover, pairs characterized by different momenta 
are incoherent (in the sense that in expectation values any product of annihilation and 
creation operators of different momenta will factorize). 

These two properties are explicit when writing the in vacuum in terms of out states 
(i.e. states with a definite out particle content). From Eq. (11), one gets (see [16, 29]) 

|0, in) = JJ<g> |0, k,m) 2 

k 

= ® ( ex P «k U<t a>™k) l°> k > out ) ® 1°' " k > out )) ■ (I 2 ) 

k 

The tilde tensorial product takes into account only half the modes. It must be introduced 
because the squeezing operator acts both on the k and the — k sectors. The definition 
of this product requires the introduction of an arbitrary wave vector to divide the modes 
into two sets. The sign of k z can be used. Notice that a rigorous definition of J| k requires 
to consider a discrete set of modes normalized with Kroneckers (that is, to normalize the 
modes in a finite conformal 3-volume). To be explicit, the two-mode state |0, k, in)2 is 
given by 

|0, k, in)2 = |0, k, in) ® |0, — k, in) , (13) 
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where |0, k, in) is the ground state of the k-th mode at the onset of inflation. The complex 
parameter z k appearing in the squeezing operator in Eq. (12) is given by the ratio of the 
Bogoliubov coefficients 

z k = % = - e - 2 ^ 1 



(1 + i2kr/ r - 2k 2 n?) 
= -l + 0(k Vr ) 3 . (14) 

The high occupation number limit corresponds to \z k \ — ► l - . 

It has to be emphasized that none of the out states in Eq. (12) carries 3-momentum. 
Hence, the distribution is homogeneous in a strong sense: at late time the 3-momentum 
operator is still annihilated by the state of Eq. (12). (This property is not satisfied by 
incoherent distributions such as thermal baths. In those cases, the 3-momentum fluctuates 
and vanishes only in the mean.) The present distribution is also isotropic since the Bogoli- 
ubov coefficients are functions of the norm k only. Finally, it is a Gaussian distribution, 
as can be seen from Eq. (12). 

To appreciate the peculiar properties of the distribution of Eq. (12), and as a prepara- 
tion for the analysis of the problem of decoherence and entropy, it is interesting to consider 
the most general homogeneous, isotropic, and Gaussian distributions of out quanta. A de- 
tailed description of these states is given in Appendix D. Their properties are completely 
specified by three real functions of the norm k (one real and one complex) through the 
following expectation values 

(0 = 0, (15a) 
(a^a<$ t )=n k 5 3 (k-k'), (15b) 
(a^ t a^ t )=c k 5 3 (k + k'). (15c) 

In the second line, n k is the mean occupation number. In the third one, the complex 
number c k characterizes the quantum coherence of the distribution. The degree of two- 
mode coherence is given by \c k \/{n k + 1/2), see Appendix C. It is bounded by 1. For a 
thermal (incoherent) distribution, one has c k = 0: no 2-mode coherence. 
In the case of pair production from vacuum, one has 

i 1 2 Zk 

nk = \Pk\ =~ : — [2> c k = a k [3 k = - ■ — . (16) 

1-Ffe| 

Therefore, when considering relevant modes in inflation, using Eq. (8) one has 



■*= m? a 10100 ' (17a) 



Cfc : 



n k + - + - 
2 n 



e l2f \ 9 k = f (k Vr f + 0(n~ 5/4 ). (17b) 



That is, the distribution which results from inflation is highly populated 2 and, more 
importantly, maximally coherent. These two properties go hand in hand for two-mode 
squeezed states. When computing the two-point function, the mean occupation number 
determines the primordial power spectrum, while the two-mode coherence of the distribu- 
tion manifests itself in the temporal coherence of the modes, as we now explain. 
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We recall how n k ~ 10 100 is obtained. For the test field £ = 4>/a, using Eq. (9), the power spectrum at 
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B. Two-point function and the neglect of the decaying mode 



After the reheating, when expressed in terms of out modes, the (Weyl ordered) two- 
point function associated with the general distribution specified by Eqs. (15) is 



G-^x^x') 



dk k 3 



sin k x — x 



k 2tt 2 &|x-x'| 



nk + 77 ) 4>k(v)<P*k(v') + Ck <t>k(ri)(t>k(v') 



c.c. 



(18) 



We have dropped the superfix 'out' on the modes for lisibility. We have chosen to work 
with the anti-commutator in order to obtain a symmetrical function in rj, rj and to analyze 
the classical limit, namely the high occupation number limit n k S> 1. 

In order to discard the contribution of the decaying modes, two conditions must be 
met. First, the sum in brackets should factorize. For a general distribution, it does not. 
However it factorizes for coherent states, see Eq. (All) in Appendix A, for distributions 
resulting from pair production from vacuum in the high occupation number limit, and, 
more generally, whenever the following inequality is satisfied 



n k + 1/2 - [cfcj 
n k + 1/2 



< 1 



Let us write 



c-k 



(19) 



(20) 



where e k and k are real. Eq. (19) is equivalent to |efc| <C 1. Then one can write the 
bracket of (18) as 



n k + 



1 



(e^faj-cc.) (V^(r/) - c.c.) + €k (j^faWMrf) + <=■<=■) ] -(21) 



One thus verifies that Eq. (19) is sufficient to factorize the two-point function. 



the onset of the radiative era is, in physical units, 

-r, \ k' A , ; . ,~, . k A h 4(krj r ) 2 k 3 h H 2 

where a" 1 — Hrj r . The power spectra of the gravitational potential \& and the gravitons hf? , are 
obtained from the above. For slow-roll inflation [24] , a mode by mode integration yields 

VA-KhG(H k \ 2 ( dine A 

Here designates the value of the Hubble parameter at horizon exit: a(rj k )H(jj k ) = k. The data from 
COBE normalize the power spectrum to be 8T k /T\ rm s — *fc/3| r ms ~ 10 -5 . The energy scale of inflation 
is thus (LpiH) 2 ~ 10 _10 e, with e = 10 -2 — 10 _1 . The mean occupation number n k must therefore obey 

n k ~ e(A r #) 4 , 

where A r = a r /fc is the physical wavelength of a mode k at reheating. Assuming that reheating happened 
at redshift ~ 10 28 , i.e. with a temperature ~ 10 14 GeV, modes corresponding to large structures today, 
e.g. lOOMpc ~ 10 21 m, have a mean occupation number ~ 10 96 . 
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The second condition for discarding the decaying mode arises from the fact that the 
mode in parenthesis, 

<j> k {r } ) = i2Im(e i6 ^T\r } )) , (22) 

will not, in general, be the growing mode. Indeed it is only for Q k <C 1 that one can 
approximate it by the sine function. 

In brief, to discard the cosines, both e k <C 1 and ^ « 1 must be satisfied. When 
considering pair creation from vacuum (8), we get e k oc l/n k and 9 k oc 1/n^ 4 . Hence it 
is perfectly legitimate to discard the cosines. In that case, the two-point function of the 
physical field £ = <j)/a reduces to 

, r £, \ 2/i is -, , > f°°dk sin(/dx — x'l) n ,, . s'm(kri) sm(kri') , . 
(0»|{«,.x),«V,x')}|0») = I T ^_ x1 " n°W lj> I, V - (23) 

The remaining statistical properties of the distribution are the power spectrum and the 
temporal coherence of the growing mode. First, the primordial power spectrum 

*?<*> - (f )' • < 24 > 

is proportional to the mean occupation number n^, see footnote 2. Here, it is scale 
invariant because the inflationary background has been approximated by de Sitter space. 
Second, the time dependent function which appears in (23) is (sin kr))/(krf), where rj is 
proportional to the scale factor given by Eq. (2b). This is how the temporal coherence of 
modes obtains from the two-mode coherence of the distribution. Consequently, the power 
spectrum V^(k;rj) at time rj, defined by the Fourier transform of the two-point function, 
is the product of the primordial power spectrum and the square of the mode function 
evaluated at that time: 

The power spectrum at a given time, as a function of k, oscillates and has zeros. This 
behaviour is a necessary condition for the existence of acoustic oscillations and anti- 
correlations in the temperature anisotropy TT and cross-correlation TE power spectra 
respectively [30]. 

Hence, once the cosine is neglected, the quantum distribution can be effectively replaced 
by a stochastic Gaussian distribution of classical fluctuations 

with locked temporal argument, and random amplitudes with variances given by 

= Uie k/ ))eff = Vl(k)5 3 (k-k'). (27) 
Being Gaussian, the effective probability distribution is simply 

*«-fimrM-Mi)- <28) 
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To avoid double counting, one must again use the tilde product which takes into account 
half the modes only, as was done in the quantum distribution of Eq. (12). This counting 
becomes crucial when computing the entropy, see Section IV. C. 

We emphasize that the reality of the field 4>(n, x) has nothing to do with the temporal 
coherence of the modes of Eq. (22). Nevertheless, once having neglected the decaying 
mode, it is true that the first equality in Eq. (27) is imposed by the reality of <f>(r), x). 
However, this should not be confused with the relation between the expectation values 
of the creation and destruction operators in the squeezed vacuum state prior to have 
discarded the decaying mode, 

(a «*t K f) = a°_t)(l + 0(k Vr f) . (29) 

The latter follows from Eqs. (15)-(17), and is the expression of the two-mode coherence 
of the in vacuum (or more generally of strongly correlated two-mode distributions, see 
Section V). The (complex) Eq. (29) guarantees that both conditions £i < 1, ^ C 1 are 
satisfied, thereby allowing to factorize the 2-point function and to discard the contribution 
of the decaying modes. 



C. Additional remarks 

First we remind the reader why a random distribution of both the sine and the cosine 
does not give rise to temporal coherence [6]. In fact such a distribution corresponds to an 
incoherent (thermal) distribution. 

Consider the incoherent distribution: 

(alpine = n k S 3 (k-k'), 

(a k a k /)^ = c k S 3 (k + k') = . (30) 

The last equation implies that temporal coherence of the modes is lost, as can be seen 
from the temporal behaviour of the bracket in the integrand of Eq. (18) : when rj = rf 
the bracket no longer exhibits oscillations in k. (This absence can also be understood by 
considering field amplitudes as classical stochastic variables rather than quantum ones. 
Writing the mode in terms of its norm and its phase [6, 31] 

k = $ k sin^r? + k ) , (31) 

one can treat $ k and # k as stochastic variables. One verifies that the distribution for 
the phase 6* k is uniform over the interval [0, 2tt]. Hence, taking the ensemble average to 
compute the power spectrum, no temporal coherence could obtain.) 

It is of value to estimate what has been neglected when discarding the cosine in Eq. 
(26). To this end, we write the field modes directly in terms of the growing and decaying 
modes 

A sm(kri) ; cosjkri) 

= + dl wr ■ (32) 
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The leading non- vanishing order in n k of their fluctuations and cross correlation are 

(9k9l)in = U k + \ - Re(c k )^j 5 3 (k - k') = 2n k (l + 0{n?j) <5 3 (k - k') , (33a) 

<44)<n = (n k + \ + Re(c k ^j <5 3 (k - k') = 0(n k 1/2 ) 5 3 (k - k') , (33b) 

({4, $£,}>i„ = /m(c fc ) <5 3 (k - k') = 0(n^ 4 ) <5 3 (k - k') . (33c) 

Notice that the cross-correlation, ({d^, g^,})i n , the Weyl ordered product (i.e. the anti- 
commutator divided by 2), is non-zero. This means that the growing and the decaying 
mode are not eigenmodes of the distribution (see Appendix C). 

When divided by the variance of g^, the variance of the cross-correlation is of order 
n fc 3 ^ 4 ~ (fer? r ) 3 ~ 10 -75 . The power of the decaying mode is even smaller. One can 
therefore safely use the effective distribution Eq. (28) in replacement of the original 
quantum distribution Eqs. (12) or (33) when calculating the power spectrum. However 
this is not the case for the entropy. 

It is illustrative to compare the in- vacuum correlations, Eqs. (33), to those of the 
incoherent distribution (30). In that case and g^ are uncorrelated Gaussian operators 
with equal variance: 

(gkgt^mc = {d k dl,)inc = + ~ k ') . 

({4,2k'}>™ = 0. (34) 
These differences with Eqs. (33) are particularly clear. 



D. Drawbacks of the simplified description 

The simplified description in terms of a statistical ensemble of sine standing waves, see 
Eq. (28), has several shortcomings. It is of value to describe them with some attention 

First, there is the discrepancy of the value of the entropy mentioned after Eq. (33). 
Second, it should be pointed out that, in the early universe, non-linear processes, however 
weak they may be, will modify the density matrix obtained by making use of free fields. 
To be able to parameterize the modifications of the power spectrum, one must return to 
two- mode distributions 3 . Indeed, the classical ensemble of sine waves cannot describe 
these modifications. It should be considered only as an effective description of the density 
matrix when the decaying mode is strictly negligible. 



3 As mentioned in the Introduction, we shall neglect non-Gaussianities when computing the entropy. In 
other words, we replace the actual (non-Gaussian) distribution by the Gaussian one which possesses 
the same moments (15) and therefore the same power spectrum. The justifications for this substitution 
are the following. First, it leads to a drastic simplification of the analysis, and all Gaussian (isotropic) 
distributions can be considered. Second, when focussing on a given wave vector, one anyway traces 
over the non-Gaussianities relating this scale to other wave vectors. Third, since non-Gaussianities are 
expected to be weak [13, 32], it is a well defined (mean field/Hartree) and useful approximation to 
replace the actual distribution by the Gaussian one associated to it. 
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The appropriate basis to describe density matrices is provided by coherent states. The 
reasons are the following. First, they constitute the quantum counterparts of classical 
configurations in phase space. Therefore they provide an adequate basis for studying the 
semi-classical limit. In particular, when non-diagonal terms are still present in this basis, 
it is the signal that some quantum coherence has been kept, see Section V and Appendix 
D.3. Secondly, they are the preferred basis in which squeezed states decohere when weak 
interactions are taken into account, see Section IV. A and [12] for further details. 



III. TWO-MODE COHERENT STATES 



When using coherent states in inflationary cosmology, one must pay attention to the 
entanglement between k and — k modes. As we shall see this leads to the notion of 'two- 
mode coherent states'. On the contrary, a naive use of coherent states assigning amplitudes 
to each mode separately will erase the entanglement and therefore suppress the temporal 
coherence of the modes. For completeness this is shown in the next subsection. 



A. Naive use of coherent states and the loss of temporal phase 

If assuming that each mode decoheres separately, the matrix density of the k, — k pair 
of modes factorizes and the modes are independent: 

p Naive = (J^-G(v)\v, k)(v, k|) (J^G(w)\w, -k)(w, -k|) , (35) 

where G(v) = n^ 1 exp ^— \v\ 2 /n^j is the normalized Gaussian distribution giving the 

appropriate power-spectrum, and where \ v, k) and \ w, — k) are coherent states constructed 
with out operators, see Appendix A for their definition. The expectation values in this 
decohered density matrix are: 

Tr(^ e 4a k ) = n k , 
Tv(p Naive a k a_ k ) = 0. (36) 

One thus verifies that Eq. (35) corresponds to the incoherent distribution of Eq. (30) 
written in the basis of one-mode coherent states. The two-mode correlations are completely 
erased and the temporal coherence of the modes is lost. In fact Eq. (35) is maximally 
decohered, since it is a tensorial product of two thermal states (because G(v) is Gaussian 
and isotropic in the complex v space). 



B. Representation of the in vacuum with two-mode coherent states 



To understand the usefulness of two-mode coherent states it is appropriate to first 
mention the following properties [33]. Consider a mode k in a coherent state \v, k). Then 
compute the one-mode reduced state obtained by projecting it on the two-mode initial 
vacuum of Eq. (13). You will get: 



(v, k|0m, k) 2 = Ak(v) \z k v*, -k) 



(37) 
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It is remarkable that the state of the — k mode which is entangled to \v , k) is also a coherent 
state. Its amplitude is given by the complex conjugate of v times z k characterizing the 
pair creation process. These facts follow from the EPR correlations in the initial vacuum 
displayed in Eq. (12). The prefactor Ak(v) is 

A k (v) = -^exp (-P-2 ) • (38) 
I a* | y 2\a k \ J 

It is the probability amplitude to find the mode k in the coherent state \v, k) given that 
we start with vacuum at the onset of inflation. 

Using the representation of the identity with coherent states, Eq. (A13), Eq. (37) 
permits to decompose the two- mode in vacuum as a single sum of two-mode coherent states. 
It should be stressed that two independent integrations over one-mode coherent states are 
necessary to describe a generic two-mode state. The fact that only one integration is 
sufficient is therefore a direct consequence of the two-mode coherence of the initial state 
|0m, k) 2 . Explicitly we have 

/j2 
-^-A k (v) \v, \s)®\z k v*, -k). (39) 

This result is exact 4 and applies even for low occupation numbers. It is the consequence 
of the coherence of the in vacuum and holds for every homogeneous pair creation process. 

Another important consequence of Eq. (37) is that the probability to find simultane- 
ously the k-mode with coherent amplitude v and its partner with amplitude w is 

V 2 , k (v,w) = \(v, k\(w, -k|0m, k) 2 | 2 = \A k (v)\ 2 x e~\ w ~ z ^ . (40) 

The second factor follows from the overlap between the reduced state in the r.h.s. Eq. (37) 
and the bra (w, — k| (we also used the coherent state overlap: |(u|-u)| 2 = exp(— \u — v\ 2 )). 

Equation (40) implies that once the amplitude of the k-mode is fixed, the conditional 
probability to find its partner in a coherent state \w) is centered around 

w(v) = z k v* . (41) 

Moreover, in the high occupation number limit we are dealing with, the spread (= 1) 
around this mean value is negligible when compared to the spread in v which is given by 
\ct k \ = n k + 1. Therefore, when computing expectation values to leading order in n k , 
the conditional probability acts as d^(w — w(v)). Hence both the real and the imaginary 
parts of w are fixed. This is how the EPR correlations in the in-vacuum translate in the 
coherent states basis: double integrations over one-mode coherent states can be replaced 
by single integrations over two-mode states. 

Furthermore, the strict correspondence between the coherent state amplitudes in the 
sectors k and — k leads, in the high-squeezing limit, to the temporal coherence of modes 



4 Notice however that the above decomposition is not unique since the coherent states are not orthogonal, 
compare with [34]. Eq. (39) has the advantage to be directly related to the detection of a quasi-classical 
configuration in the k sector. 
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since the only random variable is their common amplitude [33]. As we shall see below, 
this result is crucial as it allows to erase the quantum features of the distribution (such as 
the existence of a sub-fluctuant direction in phase space) while preserving the temporal 
coherence of modes upon horizon re-entry (which can be classically interpreted). 

To complete this analysis, and in preparation for studying decoherence, it is also in- 
teresting to explicitly write the non-diagonal matrix elements of the in vacuum density 
matrix. One has 

(v\(w\p in \v')\w'} = A 2 , k (v,w) A 2 ,k{v',w')* , (42) 
where the two-mode amplitude is 

A 2 ,k{v,w) = Ak^e-^-^l 2 e ilm ( w * z * v *K (43) 

Since the initial vacuum is a pure state and since \ak\ = ^ + 1 > 1, the above matrix 
elements do not vanish when w = z^v*, w' = z^v'* even for \v — v'\ 2 ~ n^, i.e. they do not 
vanish even for macroscopically different coherent states. (Remember that nk — 10 100 .) 

Because of this macroscopic quantum coherence, pi n does not describe a classical en- 
semble (a mixture) of coherent states. Fortunately, as we shall see below, this coherence is 
unstable to small modifications of the distribution. To describe this decoherence we shall 
proceed in two steps. First we shall consider a specific example of decohered density matrix 
and argue why it plays a critical role, namely, it separates distributions which have kept 
some quantum coherence from those which are mixtures (i.e. distributions for which the 
off-diagonal matrix elements vanish for \v — v'\ 2 > 1.) Second, we shall describe partially 
decohered density matrices in general and more technical terms. In particular we shall sort 
out the robust properties from its unstable ones. In one word, when using |/ (n^ + 1/2) 
as a parameter for characterizing the two-mode coherence of the distribution, we shall 
see that the power spectrum and the temporal coherence are robust whereas the entropy, 
the above quantum coherence, and the residual squeezing are extremely sensitive to this 
parameter. Even though the value of this parameter will be introduced from the outset, 
similar results obtain when performing a dynamical analysis. 

IV. MINIMAL DECOHERENCE SCHEME OF PRIMORDIAL SPECTRA 
A. Zurek et al. analysis and minimal decoherence scheme 

In general, it is a difficult question to determine into what mixture an initial density 
matrix will evolve when taking into account some interactions with other modes and 
tracing over these extra degrees of freedom so as to obtain an effective (reduced) density 
matrix. There exist however several cases where clear conclusions can be drawn. First, 
when one can neglect the free Hamiltonian, the preferred states (that is the basis into which 
the reduced density matrix will become diagonal) are the eigenstates of the interaction 
Hamiltonian [14, 35, 36]. This approach has been applied in [21], to primordial density 
fluctuations when the (physical) modes are almost constant because their wave length is 
much larger than the Hubble radius. The conclusion is that the preferred basis is provided 
by amplitude (position) eigenstates. However this conclusion leaves some ambiguity and 
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might lead to some difficulties. First, position eigenstates are not normalized. Second, 
and more importantly, the spread in momentum is infinite for these states. Therefore, the 
velocity field would not be well defined when the modes re-enter the horizon. Moreover, 
as pointed out in [22], some additional decoherence could be obtained as they re-enter the 
horizon. In this case, the momentum should be treated in the same footing as the position. 
To cure these problems, some finite spread in position should be introduced. One then 
needs a physical criterion to choose this spread. 

To introduce a critical spread, to identify the various regimes of decoherence, and to 
determine their range, one should use coherent states, as we now argue. (We reserve for 
a forthcoming publication a proper justification of the physical relevance of these states 
in a cosmological context.) First, as shown in [12], when considering harmonic oscilla- 
tors weakly interacting with an environment, coherent states provide the basis in which 
the density matrix decoheres. In particular, when the initial state is a squeezed state, 
there is a phase of rapid growth of the entropy which sends the system into a mixture 
of coherent states. This is then followed by a period of slower increase accompanying 
dissipation or thermalization. Second, in inflationary cosmology, fluctuation modes are 
weakly interacting harmonic oscillators [13]. Indeed, given that the relative density fluc- 
tuations have small amplitude (~ 10~ 5 ), the hypothesis of weak interactions is perfectly 
legitimate. Therefore, the modes with k' ^ k will effectively act as an environment for 
a given two-mode state. A novelty with respect to the analysis of [12] is that in cosmol- 
ogy, one deals with two-mode squeezed states. In this case one ends up with statistical 
mixtures of two-mode coherent states, but not with mixtures of one-mode coherent states. 
The procedure to define these distributions is explained in Appendix B. However, because 
coherent states are overcomplete, there is an inherent ambiguity (albeit small in the large 
n limit) in determining this distribution, c.f. footnote 4 for a related ambiguity. The 
procedure we shall adopt in this Section is to use Eq. (40) to define it. 

So unless the interactions taking place in the early cosmology are sufficiently weak 
(or sufficiently anisotropic in phase space as in the case they act only for a fraction of a 
period) so as to keep some quantum coherence, we can infer that the effective entropy (see 
footnote 3) of the final distribution should be higher than (or equal to) that stored in the 
diagonal distribution of two-mode coherent states, here after called minimally decohered. 

In the next two subsections, we shall explicitly write down this distribution, compute 
the entropy it carries, and, very importantly, show that the temporal coherence of the 
modes is kept whereas the macroscopic quantum coherence of the original distribution, 
see Eq. (42), is lost, thereby allowing a classical interpretation of the residual statistical 
properties. 



B. The minimally decohered distribution 



The minimally decohered distribution which corresponds to Eq. (42) is 




(44) 



The probability distribution Vi, k(v,w) is given in Eq. (40). 
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By having replaced the original distribution p; m = |0,m)(0, in\ by the decohered one 
Prnin-, some entropy has been introduced. Indeed, by direct computation, one gets 

f 2 

Tr(^ m i n a[o k ) = \a k \ = n k + 1 , 
Tr(p mi „a k a_k) = z k \a k \ 2 = c k . (45) 

The first line shows the occupation number increased by one unit. Thus the relative 
change is only l/n k . In the second line one finds that that the coherence term is equal 
in phase and amplitude to that of the original distribution, see Eqs. (15, 16). These two 
results imply, first, that the (relative) degree of coherence has been reduced and therefore 
some entropy has indeed been gained. Secondly, in the high occupation number limit 
we are dealing with, the two-point function of Eq. (18) is not affected by this loss of 
coherence since, for relevant modes, the relative change is of the order of l/n k ~ 1CP 100 . 
Let us emphasize that the smallness of this change follows from the fact that the phase of 
Tr(p m j n aka-k) stays equal to that of the r.h.s of Eq. (15c). Hence, the fact that Eq. (19) 
still applies leads to the preservation of the temporal coherence of modes. 

In fact, the main modification associated with the replacement of pi n = |0, in) (0,in\ by 
Pmin concerns the elimination of the non-vanishing off-diagonal matrix elements. Indeed, 
in the high squeezing limit, we have, see also Appendix D.2, 

(v\(w\p min \v')\w') 



3\n k \ 



exp-- 



A / * 



\v\ + \v I + \w\ + \w I — - I V V + w w — — V w — ^-vw 



(46) 



(We have set \z k \ = 1 at the end of calculations to simplify the expression.) The above 
result implies that along the big axis of the ellipse, see Fig. 1, i.e. for w = z k v* , the term in 
the exponential is (4/3)|f — v'\ 2 . Therefore, two coherent states separated by \v' — v\ > 1 
will no longer interfere quantum mechanically. 

Moreover, each one of the coherent states in Eq. (44) is stable under the evolution of 
the quadratic Hamiltonian, i.e., the expectation values of 0k, <9r?<Ak in these states evolve 
according to the classical equations of motion, see Eqs. (A7). Hence, the diagonal distri- 
bution (44) can safely be interpreted in classical statistical mechanics, i.e. as a (stable) 
statistical ensemble of (non-interacting) classical states. 

Finally, when computing expectation values in leading order in n k , this distribution 
can be further simplified by replacing the Gaussian factor of unit spread of Eq. (40) by 
a double delta function as discussed after Eq. (41). Using this equation, one gets the 
simplified distribution 

/d 2 v 
— \A k (v)\ 2 \v, k)(v, k| \w(v), -k){w(v), -k| . (47) 

It is now written as a single sum of two-mode coherent states, thereby making explicit 
the fact that we are dealing with a statistical superposition of two-mode states (and not 
a double superposition of one-mode coherent states) . 
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FIG. 1: Decoherence of a one-mode squeezed state in the coherent state basis. Given the 2- 
dimensional character of this sheet of paper, it is simpler to represent the phase space of a one-mode 
squeezed state than that of a two-mode squeezed state. To this end, one should first decompose 
the two-mode squeezed states we are dealing with into two copies of one-mode squeezed states, as 
exposed in Appendices B and D.4. The large red ellipse represents the 1 — a contour of a one-mode 
squeezed state in {g±, pi) space, with u> = 1. The squeezed state is characterized by the squeezing 
parameters rk et (f>k, where = sh 2 r^ and 4>k = arg(cfc)/2. The small green circles represent the 
1 — it contour of coherent states, which have no privileged direction in phase-space. When replacing 
Pin by Pmin, two modifications obtain. First, the squeezed spread in the direction of the small 
axis of the ellipse is replaced by that of the vacuum. Second, and concomitantly, the macroscopic 
quantum coherence along the big axis, see Eq. (42), is removed and replaced by that of coherent 
states. 



C. Minimal entropy and the neglect of the decaying mode 

Our aim is to compute the entropy stored in Eq. (44). The entropy of any Gaussian 
two-mode distribution, _k, can be exactly calculated [37, 38] by using the fact that 
its density matrix is unitarily equivalent to the tensorial product of two thermal density 
matrices of auxilliary oscillators, see Appendix C. These two real oscillators can be taken 
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to be the real and the imaginary parts of 0k we used in Appendix B. One can indeed write 

Pk, -k = pth, 1 ® Pth, 2 M , (48) 

where M is a unitary operator acting on the two-mode Hilbert space. The expression of 
the (von Neumann) entropy immediately follows: 

5[pk,- k ] = S[p thA }+S\p tK2 } , (49) 

where the entropy of a thermal bath with mean occupation ft is 

S [pth] = (n + 1) ln(n + 1) - ft ln(n) . (50) 

When considering distributions preserving homogeneity and isotropy, the occupation num- 
bers of the thermal matrices are equal and given by 

n k + -= ({n k + -f -\c k \ 2 \ , (51) 

where n k and c k are defined in Eq. (15). It should be noticed that the phase of c k (which 
is essential for the temporal coherence of modes) drops out from this expression. Hence 
the quantum purity is not univocally related to the temporal coherence. 

Let us apply Eqs. (50, 51) to several cases. First, for the two-mode in vacuum of Eq. 
(12), the occupation number and the coherence term are related by Eq. (16), one has 
ft k = 0, as expected. Hence the entropy vanishes. 

For the decohered matrix Eq. (44), using Eq. (45), the occupation number of the two 
thermal baths are 

n k = g (- 1 + \/8(n fe + 1) + l) ~ (52) 

where the last term is the leading order when n k 3> 1. The two- mode entropy of this 
mixture is then 

<Sk, -k [Pmin] = 2S\pth] = 2lnft k + O(l) , 
= lnn fc + 0(l) =2r* + 0(l), 

~ 1001n(10). (53) 

In the second line, we have expressed the occupation number in term of the squeezing 
parameter r k : n k = sh 2 r k . Hence, a two-mode squeezed vacuum state which decoheres 
in the two-mode coherent basis goes along with an entropy of Sk, -k = 2 r k per two-mode. 
This value is large, but not maximal. Indeed, when c k = one would have found the 
maximal value of the entropy. It is twice the above value, i.e. 

5^ k = 4-_ k = 4r fc , (54) 

or S^ c = 2r k per mode k. 

It is interesting to notice that the entropy associated with the effective distribution 
Eq. (28) of sine functions equals that of Eq. (44), up to an arbitrary constant which 
arises from the usual ambiguity of attributing an entropy to a classical distribution. (This 
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ambiguity can be lifted when introducing h to normalize the phase space integral.) Using 
this trick, the entropy associated with Eq. (28) is S e ff = 2 r k for each independent mode 
(the entropy is maximal because the state is Gaussian). However, 'for each independent 
mode' here means for each two-mode since the mode — k is no longer independent once the 
cosines have been neglected, see Eq. (27). From this equality of entropies to leading order 
in rife, we conclude that the quantum density matrix which corresponds to Eq. (28) can be 
taken to be that given by Eq. (44). A priori one might have thought that many quantum 
distributions can be associated with Eq. (28). This is not the case when imposing that 
Gaussianity is preserved and that the entropies coincide. Indeed, in the high occupation 
number limit, the entropy is an extremely sensitive function of the residual coherence. 
As we shall see in Appendix D7, this drastically restricts the space of quantum density 
matrices in correspondence with Eq. (28). 

Having obtained a lower bound for the entropy from quantum considerations, we now 
provide an upper bound for the entropy which could have resulted from processes in the 
inflationary phase. 

D. Maximal entropy and backreaction effects 

The bound follows from the fact that increasing the decoherence implies increasing 
the power of the decaying mode, see the power of the cosines in Eqs. (33) and (D17). 
One thus obtains an upper bound on the decoherence level when requiring that the power 
of the decaying mode be smaller than that of the growing mode at the onset of the 
radiation dominated era. This requirement follows from the fact that the r.m.s. value of 
the primordial fluctuations (of the Bardeen potential) cannot be much higher than the 
contribution of the growing mode from in vacuum because otherwise this would invalidate 
the whole framework of linear metric perturbations. Notice that this type of bound also 
occurs in the context of the trans-Planckian question [27, 39, 40], or whenever one is 
exploring the upper value of the deviations from the standard predictions of inflation [28]. 

Using Eqs. (32), (D17) and (D18), and expanding the sine and cosine as cos(kr/ r ) = 
1 +0(n~ 3 / 2 ) and sm(krj r ) = 0(n -1 / 4 ) where we used oc (fcr/ r ) -4 S> 1, the condition on 
the power of the decaying mode stated above reads 5k < C^n^ 1 / 2 ), or equivalently 

n k <n k /A , S k ,_ k <3r fc + 0(1). (55) 

If no further decoherence is added in the radiation dominated era, this is the maximum 
amount of entropy stored in the primordial spectrum. Notice that when evaluated at 
recombination, the two-point function is still unaffected by this modification of the co- 
herence because at that time the decaying mode has still further decreased. Indeed the 
relative modifications are then of the order of n^ 1 ^ 2 ~ 10~ 50 . 

From this analysis we (re-)learn that the dynamical processes occurring in the inflation- 
ary era which are compatible with a linearized treatment of perturbations will, in general, 
leave no significant imprint on the spectra at recombination. 
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V. HOMOGENEOUS GAUSSIAN DISTRIBUTIONS WITH HIGH 
OCCUPATION NUMBERS 

The purpose of this section is to identify the quantum and classical statistical properties 
of homogeneous Gaussian distributions when the two-mode coherence is not perfect. In 
particular, we shall determine the two interesting ranges of decoherence level. First there 
is a narrow range wherein modifications of the initial density matrix are so small than the 
distribution keeps quantum properties. Second there is a much wider range wherein the 
distribution, on the one hand, has lost its quantum properties so as to be interpretable in 
classical terms, but, on the other hand, has preserved the temporal coherence of modes. 
We shall only present the results. The interested reader will find the technical details in 
Appendix D. 

The most general statistically homogeneous and isotropic Gaussian distribution is char- 
acterized by its second moments given in Eqs. (15). They are characterized by three real 
quantities. However when focussing on the coherence of the distribution, only one param- 
eter matters, namely that governing the norm of the coherence term. In this paper, we 
shall not discuss the possible modifications of the spectrum which are due to a change of 
the phase of c&. Firstly because this change induces no entropy change. Secondly because 
these modifications have been discussed in the context of the trans-Planckian question 
[27, 39, 40]. In that case indeed the phase shifts result from a change of the Bogoliubov 
coefficients which preserve the purity of the state. Therefore those changes are orthogonal 
to the ones considered in the present paper. 

A convenient parametrization of the coherence level is provided by 5k which is defined 

by 

\c k \ 2 = (n k + l)(n k - 5 k ) ■ (56) 

5k ranges between and n^- They are three characteristic values: 5k = obviously 
corresponds to the pure squeezed state: the original in-vacuum, 5k = 1 corresponds to the 
critical value above which density matrices loose their quantum properties, see Appendices 
D3-6, and 5k = n k corresponds to the thermal case with maximal entropy. 

A. Quantum and classical properties 

The analysis of Appendix D reveals that the statistical properties of the distribution 
fall in two categories, according to whether they are robust or fragile when the distribution 
is slightly perturbed. By definition the "classical" properties of the distribution, i.e. the 
physical properties observable in the classical regime, are the robust ones. 

1. Classical properties 

By robust, we designate the statistical properties of the original distribution pi n which 
are preserved in a wide class of partially decohered distributions. More precisely, to 
leading order in rik, these properties are unaffected by a modification of the coherence 
term A|cfc|/|cfc| <C 1. Thus one has: 

classical robust unaffected when 5k <C rik ■ (57) 
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Since the phase of Ck is unchanged, 5k <C guarantees that the power spectrum and the 
temporal coherence of the modes are robust properties of the state. As we shall see below, 
(together with the Gaussianity which is here exactly preserved) they are the only ones. 

For high occupation numbers, the condition 5k < % is a very mild constraint. Indeed, 
robust properties will be lost only when the state is close to a thermal state (and thus 
completely incoherent). When starting from a squeezed state, this regime can only be 
reached dynamically in the presence of strong non-linearities. Given the amplitudes of 
primordial fluctuations ~ 1CP 5 , this possibility seems excluded in inflation. 

It is interesting to notice that the condition 5k <C nk can be interpreted in two comple- 
mentary ways. First, with the mean values. In this case, 5k <C n& gives (|cfc| —nk)/nk <C 1. 
This coincides with Eq. (19) which is the condition for the two-point function to factorize. 
Alternately, one can consider the distribution itself (see Appendix D.2). One gets 



(v, k|(u>, — k|/3|tt>, — k)|u, k) oc exp 



w — w 



rik + l l + S k 



(58) 



where w = CkV*/(n + 1) designates the most probable value of w given v. It generalizes 
Eq. (41). When 5k <C rik the spread of w around w is negligible with respect to the 
spread in v, the power (= rik + 1). Thus the condition 5k <C nk means that the modes k 
and — k are still tightly correlated in phase and amplitude, as they were in the original 
distribution, see discussion bellow Eq. (40). These distributions thus are (to leading order 
in n) statistical mixtures of two-mode coherent states \v)\w) as in Eq. (47). 



2. Quantum properties 

As discussed in Appendix D, the "quantum features" of the distribution are the fol- 
lowing: 

• The correlations between macroscopically separated semi-classical states, see dis- 
cussion after Eq. (42) and Appendix D.3, 

• the existence of a sub-fluctuant mode in phase-space, see Appendix D.4, 

• the non-P-representability of the state, see Appendix D.5, 

• the non-separability of the state and the violation of Bell inequalities, see Appendix 
D.6, and 

• an entropy smaller than that of the minimal scheme, see Eq. (53) and Appendix 
D.7. 

The important fact is that all these properties are erased when 5^ > 1, independently 
of nk when nk 3> 1. One can thus say 

quantum <=^ fragile <=^ erased iff 5k > 1 . (59) 

Let us emphasize the important lessons. The first is that these quantum properties 
are still present when 5k < 1. Therefore 5k = 1 is indeed the critical value from which 
quantum properties are lost. 
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Second, the density matrix which corresponds to 6k = 1 is that given in Eq. (47). From 
Eq. (D23), one clearly sees that it is the least decohered separable distribution. 

Third, for n& 3> 1, the quantum properties are very sensitive to a slight perturbation 
of the system since 5k = 0(1) corresponds to a relative change of order l/n^ of the power 
spectrum Eq. (21). Therefore, there is a wide range of 5k (from 1 to nj, with 7 < 1) 
where the robust properties are still unaffected (to leading order in rik) while the quantum 
properties are completely erased. 

Fourth, for 5k < 1, the product of the residual squeezed spread times the residual 
width of the correlations between semi-classical states in the orthogonal direction stays 
approximately constant whereas it increases linearly with 8k for 5k 3> 1, see Eq. (D16) 
and Fig 3, see also [41]. This means that these two spreads are equivalent expressions of 
the residual quantum properties of the state. 

Fifth, the entropy (which is a measure of the degree of purity of the state) is also very 
sensitive to the level of coherence, see Fig. 5. Therefore it can be used to parameterize 
the level of coherence and to order different distributions, see Appendix D7. 

B. Discussion 

First, the above results show that the following statements are equivalent in the high 
occupation number limit, for decoherence levels obeying 5k <C n& : 

• The distribution of the modes k and — k is strongly correlated in phase and am- 
plitude, where strong means that, given a value v of the amplitude of the mode k, 
the conditional probability to find — k with amplitude w is centered around a value 
w(v) with a spread = 1 + 5k <C n^, see Eq. (58). 

• All realizations of the modes contributing to the two-point function have the same 
temporal phase. This phase is fixed by arg(u)) of Eq. (58) and is not affected by 
the level of decoherence. Indeed the modifications of arg(u)) preserve the quantum 
purity of the distribution and are orthogonal to those considered here. The absence 
of phase shift guarantees that the decaying mode stays subdominant for 5k <C n^, 
in the sense that its power (= 5k/2 + 0(l/n)) stays <C rifc, see Eq. (D17). 

• The two-point function factorizes in a product of two growing modes. 

Second, both the power spectrum and the temporal coherence are fairly robust, in the 
sense that they are lost only for 5k — n-k, i.e. through a significant regeneration of the 
decaying mode. They are the statistical properties of the state which can be interpreted 
classically. The fact that the criterion of robustness singles out certain statistical properties 
should be conceived as a definition of the classical (statistical) information contained in 
the quantum density matrix, in agreement with Gottfried's analysis [14]. 

In the case where 8k < 1, the distribution cannot be interpreted in classical terms. 
Nevertheless, when probing (or observing) the momenta of this distribution with a relative 
uncertainty larger than 1/wfc, the above robustness prevents the existence of signals telling 
us that some quantum properties are still present. Indeed, quantum properties of Gaussian 
distributions can onlybe put into evidence in two cases. Either one must be able to measure 
the momenta of distribution (15) with a precision higher than l/n/<;. This amounts to be 



23 



able to measure the amplitude of the sub-fluctuant mode and hence that of the decaying 
mode, see Eq. (D17). Or one must have a set of non-quadratic operators. A first exemple 
of such operator giving rise to violations of Bell's inequalities can be found in [42], see 
also Appendix D6; a second is discussed in [43]. The consistency of these two alternatives 
is understood when noticing that the measures involved in the second case amount to 
distinguish configurations with a precision higher than l/n&. 

These considerations imply that the non-linear processes giving rise to structure for- 
mation will also be insensitive to the existence of a residual degree of quantum coherence. 
(In this we disagree with the claim [44] according to which only the decohered part of 
the power spectrum will participate to structure formation). The question of the residual 
quantum coherence (i.e. the value of 5k if 5k < 1) seems therefore of purely academic 
interest. 

VI. CONCLUSION 

In this paper we analyzed partially decohered distributions susceptible to describe the 
primordial fluctuations in inflationary scenarios. 

First we have identified the parameter governing the level of decoherence and computed 
the associated entropy, the residual value of the squeezing, and the residual coherence 
encoded in off-diagonal matrix elements. 

Second, we showed that the correlations of modes in the initial distribution and its 
decohered versions, which is at the origin of the temporal coherence of the mode at hori- 
zon re-entry, is properly taken into account by considering superpositions of 'two-mode 
coherent states'. 

Third, we discussed how, in the large occupation number limit, the fragility of quantum 
properties of the distribution to small modifications lead to a critical decoherence scheme 
above which all quantum properties are erased. The associated entropy is very large (and 
equal to 1/2 the thermal entropy) even though the relative change of the power spectrum 
is only l/n^. For 8^ > 1, the remaining statistical properties can be interpreted classically. 

The open question concerns the calculation of the residual degree of coherence (the 
value of 5k)- This is a difficult question which calls for further work. During the infla- 
tionnary period, the dynamics is indeed not trivial because of the conjunction of mode 
amplification and non- linear effects. 

Acknowledgments: We are grateful to Serge Massar for having explained to us the 
meaning of decoherence by drawing Fig. 1 many years ago in Brussels. We are also grateful 
to Dani Arteaga, Claus Kiefer, Jihad Mourad and Alexei Starobinsky for useful remarks. 

APPENDIX A: COHERENT STATES 

This appendix aims to present the properties which we shall use in the body of the 
manuscript. For more details, we refer to [45-47]. 

Coherent states (of a real oscillator) can be defined as eigenstates of the annihilation 
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operator: 

a\v) = v\v) , (Al) 
where v is a complex number. In Fock basis it is written as 

\v)=e-—Y,^\n), (A2) 



n=0 v 71 

where the exponential prefactor guarantees that the state is normalized to unity {v\v) = 1. 
They are also obtained by the action of a displacement operator on the vacuum : 

\v) =D(v)\0) =e v * & - v&t \0) . (A3) 

The first interesting property of coherent states is that they correspond to states with 
a well defined complex amplitude v. Indeed, by definition (Al), the expectation values of 
the annihilation and creation operators are 

(v\a\v) = v , {v\a f \v) = v*. (A4) 

Thus the mean occupation number is 

(v\aJa\v) = \v\ 2 . (A5) 

It is to be also stressed that the variances vanish: 

Aa 2 = (v\a 2 \v) - (v\a\v) 2 = 0, Aa f2 = {v |a f2 |t>) - (v\a^\v) 2 = . (A6) 

From these properties one sees that the expectation values of the position and momen- 
tum operators (in the Heisenberg picture) 



are 



q(t) = {v\m\v) = \J^(ve-^ t + v*e i ^) = ] J^\v\cos(u;t- ( t> v ), 



p(t) = (v\p(t)\v) 



i\J^-{ve~ iujt - v*e lujt ) = -VmJ\v\ sm(cot - <f> v ) = 8 t q{t) .(A7) 



We have used the polar decomposition v = \v\e l( ^ v . These expectation values have a 
well defined amplitude and phase and follow a classical trajectory of the oscillator. This 
is due to the "stability" of coherent states under the evolution of the free Hamiltonian 
2H = p 2 + uj 2 q 2 : if the state is \v) at time to, one immediately gets from (A2) that at a 
later time t, the state is given by \v(t)) = \ve~ lu} ( t ~ t °)). Notice that the variances of the 
position and the momentum are 

Af = -, Af = -. (A8) 



25 



They minimize the Heisenberg uncertainty relations and are time-independent. Hence, in 
the phase space (q,p), a coherent state can be considered as a unit quantum cell 2-kH in 
physical units (see also (A14) for the measure of integration over phase space) centered on 
the classical position and momentum of the harmonic oscillator (q(t),p(t)). In the large 
occupation number limit \v\ S> 1, coherent states can therefore be interpreted as classical 
states since Aq/^q 2 + p 2 /to 2 = Ap/^to 2 q 2 + p 2 = l/2\v\. This is a special application 
of the fact that coherent states can in general be used to define the classical limit of a 
quantum theory, see [47] and references therein. 

One advantage of coherent states [45] is that the calculations of Green functions re- 
sembles closely to those of the corresponding classical theory (i.e. treating the fields not 
as operators but as c-numbers) provided either one uses normal ordering, or one considers 
only the dominant contribution when \ v\ 3> 1. We compute the Wightman function in the 
coherent state \v) 

G v (tX) = (v\q(t)q(t')\v) 

= (: mW) :>„ + 1;^-'') , (A9) 

where we have isolated the contribution of the vacuum. The normal ordered correlator is 
order \v\ 2 : 

(■ q{tW) -)v = -Re \{a 2 ) v e~ w ^ + {a)a) v e w ^ 

UJ L 

= — \v\ 2 cos(u)t - <f> v ) cos(ujt' - <f) v ) = q(t)q(t') . (A10) 

We see that the perfect coherence of the state, namely |(aa)„| = {a)a) v is necessary to 
combine the contributions of the diagonal and the interfering term so as to bring the 
time-dependent classical position q(t) in Eq. (A10). 

The wave-function of a coherent state in the coordinate representation is given by 

where v = (ujq + ip) / V 2ujH. This follows from the definition (q\a\v) = v(q\v). From this 
equation one notes that two coherent states are not orthogonal. The overlap between two 
coherent states is 

(v\w) = exp ^v*w — t^M 2 — ^l^l 2 ^) ' (A12) 

Nevertheless they form an (over)complete basis of the Hilbert space in that the identity 
operator in the coherent state representation {\v}} reads 

1 



h-^\v){v\. (A13) 



The measure is 

d 2 v d(Rev)d(lmv) dqdp 



(A14) 



7r 7r 2irh 

The representation of identity can be established by calculating the matrix elements of 
both sides of the equality in the coordinate representation {\q}}, with the help of (All). 
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APPENDIX B: APPLICATION OF ZUREK & AL. RESULTS TO THE 

COSMOLOGICAL CASE 

In this appendix we show that Eq. (44) is the minimally decohered distribution in the 
sense of [12]. To this end, we decompose the complex mode 4>k into two real oscillators 4>\ 
and 4>2 given by its real and imaginary parts. Since the two-mode Hamiltonian is quadratic 
and Hermitian, it splits into the sum of two identical one-mode oscillator Hamiltonians for 
01 and 4*2 separately. Notice that the annihilation operators of these two real oscillators, 

ai = (a k + a_ k ) /\/2, a 2 = -i (a k - a_ k ) /\/2, (Bl) 

mix k and — k annihilation operators. Hence they can easily enforce correlations between 
these two modes. However, for our decoherence procedure to be valid, as noticed in [21], 
it is important that the interactions do not mix <j}\ and 4>2, so that they do not break the 
homogeneity. 

A two-mode squeezed state |0m, k) 2 can always be written as the tensorial product 
of the two one-mode squeezed states [48]. In our case, the one- mode squeezed states are 
those of the oscillators 1 and 2 because the Hamiltonian separates. Thus we have 

|0m,k) 2 = |0m,l) <g> |0m,2) . (B2) 

The one-mode squeezed states are governed by the same parameter z/2 : 

|0m,l) = ^£(f) n ^:|2n,l>. (B3) 

The same expression holds for the ket |0in, 2). 

The overlap of this one-mode squeezed state with a one-mode coherent state is 

(v, l|0m, 1) = -= exp (-\vi\ 2 /2 + zv\ 2 j2 \ . (B4) 

i / 1 rv I V / 



According to [12], when taking small interactions into account, the density matrix of a 
one-mode squeezed state will preferably decohere into the mixture 

—^Viivi) \vi)(vi\ . (B5) 

where the statistical weight is given by the probability to find a coherent state starting 
from the initial state: 

Vi{vi) = \(v, l|0m,l)| 2 = j^exp (~M 2 + i?e(z^ 2 )) . (B6) 

It should be noticed that this choice of the probability is not unique. The procedure 
adopted here consists in taking the so-called Q-representation of the original density ma- 
trix (see Appendix D.l for its definition) and to treat it as a P-representation to define the 
the decohered matrix density. A slightly more coherent distribution is given at the end 
of Appendix D5. Since the ambiguity results from the fact that coherent states are over- 
complete, no differences show up to leading order in n^. Hence, the various distributions 
belong to the same class: they have the same entropy, see Appendix D7. 
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To complete the proof we need to show that product of two one-mode coherent states 
1 and 2 is also the product of a coherent state for the k and — k modes. This is the case: 



\vi)®\v 2 ) = D 1 (v 1 )D 2 (v 2 )\0in,l) <g> |0m,2) , 

= £> k (u)£>_ k («;)|0OT,k) <g> |0m,-k) , 
= \v, k) <g) -k) , 



(B7) 



where the displacement operator D is defined at Eq. (A3), and where the amplitudes are 
related by 



v = 



vi + iv 2 



w = 



V\ — IV 2 



V2 



(B8) 



Finally, the product of the probabilities (B6) gives the probability Eq. (40). Performing 
the change of variables of integration from (v±, v 2 ) to (v, w) completes the proof. 

APPENDIX C: THE COVARIANCE MATRIX 

Let us consider a given pair of modes (k, — k) in a Gaussian state (15). We prove the 
relations (48-51). This is best done with the covariance matrix of the state. We also define 
the sub- and super-fluctuant modes. 

1. Definition of the covariance matrix 



To conform ourselves to the usual definition of the covariance matrix [49] , we introduce 
canonically conjugated and hermitian 'position' and 'momentum' variables for each mode, 
i.e. d k = (<? k + ip\) / y/2, and similarly a_ k = (<?_ k + ip-k) / \[2. All the information of the 
Gaussian state is encoded in the covariance matrix. In terms of the 4 vector 



/ <?k \ 
Pk 
9-k 
VP-k/ 



(CI) 



the covariance matrix has the matrix elements Cjj = Tr (^p | Q, Cj}) > where { , } is the 
Weyl ordered product, i.e. the anticommutator divided by 2. For the states with variances 
(15), one has 



( n+\ 



C 



\ 







n+ ■ 

Ci 



Ci 

n + 



Ci 



n + \) 



(C2) 



where c r = Re(c) = —\c\ cos 29, and Cj = Im(c) 



\c\ sin 26, see Eq. (20). This covariance 



matrix is not bloc diagonal since the modes k and — k are correlated. As we shall see in 
Appendix D6, this will lead to the notion of non-separability. 
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2. The sub- and super-fluctuant modes 

It is appropriate to return to complex modes to express the eigenmodes and the eigen- 
values of the covariance matrix (C2). They are given by the rotated modes: 



$ k = —(a k e- i9 + al k e +ie ), 



n k 



(C3) 



Their expectation values are 



<{*k,4»=(n + ^)+M, <{n k ,njj> = (n + £)-|c|, <{<& k ,n£}> = o. 

(C4) 

The first two r.h.s. expressions give the eigenvalues of (C2). They are twice degenerate 
because of homogeneity. In the case of pair creation from the vacuum, n = sh 2 (r) and 
|c| = chrshr, they reduce respectively to e 2r /2 and e~ 2r /2. They are therefore called 
the super-fluctuant and sub-fluctuant modes respectively. It is worth noticing that their 
product gives 1/4 as in the vacuum. This is the expression of the purity of the state. 

3. Proof of Eq. (48) 

An efficient way to prove Eq. (48) and Eq. (51) consists in using the the real and the 
imaginary parts of k , see Eq. (Bl). The passage from (k, — k) (1, 2) is performed 
by making a global rotation R g of angle 7r/4. By global we mean a 4 x 4 transformation 
mixing k and — k sectors. Explicitly, one has 



(C5) 







( 


1 





1 


°\ 


/ & \ 




1 





1 





1 




TT 2 /Vk 






-1 





1 





9-k 


\ -Vk(j>2 ) 




\ 





-1 





1/ 





The covariance matrix becomes bloc diagonal: 



C< 



C(l,2) - RgCR g 



\ 



2 H - 


Ci 








Ci 
















T^k ~i~ 2 ^ r 


Ci 













■ (C6) 



The state is a tensor product p = p± <£> p2 and the matrices p\ and p2 coincide. These 
properties follow from the homogeneity of the state. 

The transformation Eq. (48) amounts to bring C under the form 



C = MTM T , T 



n+\)l. 



(C7) 



where the matrix T is the covariance matrix of the product of the two thermal density 
matrices pth,i ® Pth,2 m Eq. (48). The matrix M is the product of three special matrices 



M = Si oc (k) Rloc(<P) Rg ■ 



(C8) 
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First, the global rotation R g of (C5). Second, the product of 2 x 2, i.e. local, rotations, 
Rloc(<p) = R(ip) ® R(f), where 

r-,/ \ f cos 09 sin a? \ 1 . . 7T 

= ■ > V = 7T ar g (c) = 77 + 9, C9 

V - sin (p cos a? y r 2 w 2 v ; 

brings C(i,2) into a diagonal form: diag(A + , A_, A_, A + ) with A± = n + l/2±|c|, see (C4). 
Third two local squeezing S(k) <8> S(—k) 

S(k) = ( ^ g _° K/2 ) , th« = -\c\/(n + 1/2) , (CIO) 

rescale the eigenvalues to a common value n. The latter is fixed by the conservation of 
the determinant: 

[n + £) 4 = det C = [in + I) 2 - |c| 2 ) ' . (Cll) 

APPENDIX D: GENERAL DESCRIPTION OF HOMOGENEOUS GAUSSIAN 

STATES 

We review the statistical properties of two-mode (or bi-party) Gaussian states. They 
have received much attention because of their importance in the contexts of Quantum 
Optics and Quantum Information, see for instance [50] for a comprehensive review. 

Since in cosmology one deals with statistically homogeneous and isotropic distributions, 
one needs to consider only a sub-class of Gaussian states. These symmetries considerably 
simplify the discussion since the density matrices are characterized by only three real pa- 
rameters given in Eqs. (15). It should be noticed that the squeezed states (12) (including 
the case of zero squeezing, i.e. the vacuum) are the only pure states. They are char- 
acterized by only two independent parameters. This is because squeezed vacuum states 
minimize the Heisenberg uncertainty relations, see Eqs. (D5) and (D6). The thermal state 
is the most decohered distribution and corresponds to c = 0. As in the body of the text, 
we shall parameterize the level of decoherence by 5 defined by 

|c| 2 = (n + l)(n-<5), (Dl) 

and leave the phase of c unchanged. Nevertheless, for completeness, we consider in Ap- 
pendix D.8 the effects of a smearing of the phase of c. 

1. The Q-representation of a density matrix 

To proceed, it is convenient to work with a representation of the density matrix p in 
phase-space. We choose the so-called Q-representation Q p . We recall that it is well defined 
for every hermitian, positive definite operator (density matrix). For a two-mode state, Q p 
is the expectation value of p in a pair of coherent states 

Q p (v,w) = {v, k\(w, -k|/3 k> _ k |u>, -k)\v, k) . (D2) 
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It is a function over phase space with the convention v = (k^ + iir^) / V 2fc , and w = 
(fc^_k+ivr_k)/ V%k, see Appendix A. It is remarkable that although (D2) is an expectation 
value, the knowledge of the function Q p (v,w) is completely equivalent to that of p. In 
particular the information about the off-diagonal matrix elements (v\(w\p\w')\v'), v ^ 
v', w / w' is encoded in (D2) [10, 51]. This is due to the redundant (overcomplete) 
character of the basis of coherent states 5 , see Appendix A. 

In addition, for Gaussian states, the knowledge of the expectation values (15) is all 
one needs to write the Q-representation of the density matrix. The general method of 
characteristic functions is given in every textbook of Quantum Optics, see for instance 
[10, 50, 51]. The case of homogeneous distributions can be worked out by hand, and one 
finds 



Q p (v,w) = ^-exp 



-— ((n + 1)(H 2 + H 2 ) - 2Re(c*vwfj 



(D3) 



where 



A = (n + l) 2 -|c| 2 = (n + l)(l + 5). (D4) 



One verifies that (D3) gives back (40) when considering two-mode squeezed vacuum states, 
i.e. A = n + 1. 

A necessary and sufficient condition for (D3) to be the expectation value of a density 
matrix (a positive hermitian operator) is 



n{n+ 1) - \c\ z > 0. (D5) 

This is easily seen from the fact that p is unitarily equivalent to a tensor product of thermal 
states, see Eq. (48) and Appendix C. From (51), the positivity condition n > gives (D5). 
It should be stressed that (D5) is the Heisenberg's uncertainty relations 6 , 

(akaj^a^a-k) > (a k a_ k )(a| c al k ) . (D6) 

Notice that the lower bound is reached for the squeezed states (12). 

Incidentally, the positivity condition (D5) furnishes a clear way of distinguishing clas- 
sical phase-space distributions (i.e. Probability Distributions Functions, the states being 



5 Let us work with one-mode states for simplicity. Consider the expansion in the Fock representation of 
a (bounded) operator O = Yl \n}Onm{m\. Using (A2), one sees that the function 

n,m— 

0(v* , v') = exp 

is an entire analytic function of both variables (v*,v'). A theorem for functions of several complex 
variables states that if 0(v*,v') vanishes on the line v' = v* , then it vanishes identically [51]. Consider 
then two density matrices pi and p2 with identical diagonal matrix elements in the basis of coherent 
states. Setting, O = p\ — f>2, and using the previous theorem, one concludes that p\ — p^- 

6 We recall that the general form of the uncertainty relations is, for all hermitian operators A and B, 

AA 2 AB 2 > /4. To get (D6) use A = k 0+ and B = 7r k 7r£. 



2 [m 



+ \vf) (v\6\v')= jr o n 



<yr<yr 

Vn!m! 
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points in phase-phase), and phase-space representations of quantum states. Were Q p (v, w) 
be a PDF, the positivity criterion would have been A > which is less restrictive than 
(D5). 

The off-diagonal matrix elements are [10, 50] 



(v\ (w\p\w'} \v' 



A" exp 

x exp 



1 

~2 

1 - 



i2.ii/2, ,2,, ,/2^ 

v\ + \v\ + \w\ + \w\ 



n + l 
A 



/ * / , * A i ^ / * / / i * *\ 

[y v - J rww)-\- — [cvw + cv w ) 



,(D7) 



and the density matrix itself is 
1 



P = ^exp 



: exp 



n+l / + .. 





"c* A A " 


: exp 





where ": : " stands for normal ordering. This completes the proof that the knowledge of c 
and n allows to reconstruct the quantum distribution. 



2. Factorisability of the two-point function, the classical statistical properties of 
the state, and two-mode coherent states 



The function Q p can be re- written as the product of a marginal probability Q v {v) = 



f—^- Q(v,w) and a conditional probability Q w \ v (w) 
Q p (v,w) = Q v {v) x Q w \ v (w), 



n + l 



exp 



(n+l) 



1 + 5 



exp 



\w — w(v)\' 
1 + S 



(D8) 



Both Q v and Q w \ v are normalized. The quantity w(v) = cv* /{n + 1) is the conditional 
mean value of w given v. One recovers Eq. (41) for a squeezed vacuum state. The 
conditional width around this conditional mean value is 7 



A v w 2 



w 



\w(v)\ z 



1 + 5. 



(D9) 



The means and variances can either be read from (D8) , or calculated using v and w as random variable 
whose PDF is (D8). One immediately gets 

„ .. fd 2 vd 2 w _ . . „ ,. . l2 u fd 2 vd 2 w, l2 _ , . 

«"»= / uQ P (^,w)=0, «k| 2 ))= / \v\ 2 Q p (v,w) =n + \. 

J IT TT J IT TT 

Similar equations hold for the averages (( w )) =0 and (( \ w\ 2 }} = n + 1 since (D3) is symmetric in v and 
w. The conditional means and variance of w are calculated with the conditional PDF, 

(w) = J—^- wQw\v{w) = u>(v) , A v w 2 = (|™| 2 ) - |w| 2 = y-^— (IH 2 - l™| 2 ) Qu>|u(w) = 1 + 5 . 

Notice that this spread is independent of the value of v. This guarantees the temporal coherence in the 
mean Aw 2 = (( |w| 2 }} — (( \w\ 2 )) = 1 + 8, see discussion below (Dll). For subtleties concerning operator 
ordering using phase-space representations of the density matrix, we refer to [10, 51]. 
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In the case of two- mode squeezed states (5 = 0), this is exactly 1. (Thus the second term 

2 I — ~| 2 

of the exponent is the overlap of coherent states = e~' w w < , and one recovers 

again the pure state Eqs. (39, 40)). For the decohered distribution (44) with parameter 
5 = 2, one has Aw 2 = 3. This value should be compared with the corresponding width of 
the centered distribution in v, that is Av 2 = ((\v\ 2 )) = n + 1. The ratio of the spreads is 

Al " 2 1 + S (D10) 



Therefore we obtain 



A...2 



,W 



Av 2 



Av 2 n + 1 



<1 <5<n. (Dll) 



First, when it is verified, the conditional probability to find w given v is sharply peaked 
around the conditional value w. Therefore the modes k and — k are strongly correlated in 
phase and amplitude as discussed below (40). Consequently, the power spectrum and the 
temporal coherence, see Eq. (25), are preserved whenever decoherence processes induce 
5 <C n. Put differently, whenever (Dll) applies, these statistical properties of the distri- 
bution are as well described by statistical mixtures of entangled two-mode coherent states 
\v) \w), as in Eq. (47). 

Second, (Dll) is the condition of factorizability of the two-point function Eq. (19). 
What we shall learn is that this is a rather mild constraint which includes a wide range 
of density matrices. Indeed as we shall see below, this condition is satisfied by density 
matrices which have completely lost their quantum properties. 



3. Macroscopic correlations 

As discussed after Eq. (42), the pure squeezed vacuum state cannot be interpreted as a 
classical distribution because macroscopically separated (n > \v — v'\ 2 3> 1) semi-classical 
configurations are correlated. We now show that this property is very sensitive to the 
value of 5. 

The width of these correlations is governed by the coefficient of the cross terms 
v*v', w*w' in Eq. (D7), i.e. 1 — (n + 1)/A = 5/(5 + 1). It belongs to the interval 
[0, n/(n + 1)]. The lower and upper bounds are reached for squeezed vacuum states and 
thermal states respectively. To better understand its physical meaning, it is instructive to 
pose and examine these two limiting cases. The former is presented in Section III.B. For 
a one-mode thermal state, one has 



(v\p\v')\ = ^pyexp 



f -i , \ \ /|2 1 , , ,,2 

"7 H rr \\v-V\ -— — —r\v + v'\ 

4 V n + l J 4(n + 1) 



The first term in the exponent is bounded by 1/2, which means that the two different 
semi-classical configurations are uncorrelated. The second term is the power spectrum. 

We are now ready to interpret the results for (D7). To simplify the discussion, we 
calculate the width of the off-diagonal matrix elements (D7) when w = w = cv* j(n + 1) 
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and w' 



w 



cv'*/(n + 1), i.e. for the most probable value of the conditional amplitudes 



of the — k mode. Then one has 



\(v\(w\p\v')\w'}\ = — exp 



+ 



tf(5) 2(n + l) 1 + 5 



I / 1 2 I , n2 

"2^)~~4faTl) 

1 + 5 
n + 1 



(D12) 



The width cr,, on the line u — v' is written as a function of <5. Its inverse is plotted on Fig. 
2 for a value of n = 100. For large occupation numbers, (7^(5) drops from its maximum 
value 2(n + 1) (for the two-mode squeezed state) to 1 over a range of 5 = 1 + 0(l/n). In 
conclusion, 



uncorrelated semiclassical states 



a v < 1 



5 > 1. 



(D13) 





FIG. 2: XTie inverse width a~ 2 of the off-diagonal matrix elements as a function of the degree of 
coherence of the distribution for n = 10 1 (lower, blue), n = 10 2 (middle, red), and n = 10 3 (upper, 
black). Left panel, as a function of C — S/n, one clearly sees the extreme sensitivity of a v with 
respect to 5 for small values. On the right panel, the same spread as function of x = arctan(ln((5)). 
This further shows that, for all values of n 3> 1, the critical regime is centered around 8=1, i.e. 
the 'grain' of coherent states. 



4. Sub-fluctuant mode and the power of the decaying mode 

From the second equation in (C4), we see that the sub-fluctuant mode behaves as 

'5 , (1 + 5) 2 



<{n k , njj> = [- + 



8n 



+ 0{n 



(D14) 
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Therefore one looses the squeezed spread, i.e. ({Ilk, n^}) > 1/2, when <5 > 1. In conclu- 
sion, one has 

no subfluctuant variable 5 > 1 . (D15) 

It is interesting to notice that the value of the quantum coherence a v and that of the 
squeezed spread are correlated. Indeed, their product obeys 

°l <{n k , 4 }) = ^ (i + o{± S -fj , (D16) 

for < 5 <C n. (This result is non-trivial for 5 < 1/n.) In Fig. 3, we have represented the 
plot of this product for the whole range of 5. It is growing with 5, the lower bound 1/4 
being reached for pure (squeezed) states. More importantly it stays almost constant for 
5 < 1 before growing linearly with 8 for 5 S> 1. 

The first regime corresponds to the loss of the quantum coherence whereas the second 
corresponds to the linear growth of the power of the subfluctuant mode (since stays 
~ 1). At this point it is important to notice that this second regime also corresponds to a 
linear growth of the power of the decaying mode. Indeed, the powers of the subfluctuant 
mode and the decaying mode are related by 

<{dk, 4» = <{ n k> n t» + M(l-cos20) 

= ({II k , IT*}) + 0(n-V2) = * + 0(n -i/2) , (D17) 

where 6 is defined in Eq. (20). In the second line, we have used the inflationary value of 
8 (oc n~ 3 / 4 ) and the fact that 6 is unchanged when increasing 5. 
Similarly, the cross-correlation is 

<K, 2 k }> = " sm(29)\c\ oc n 1 / 4 + + 0(n" 2 )) . (D18) 

One sees that the cross term is unaffected as long as 5 <C n. 



5. P-representability of partially decohered density matrices 

We now discuss the notion of P-representability of a density matrix [46]. It is an intrin- 
sic property of states which enlightens the distinct behaviours of quantum and classical 
distributions. 

A one-mode density matrix is said to be P-representable if it can be written as a 
mixture of coherent states 

p P . repr = f— P{v) \v, k) (v, k| , (D19) 



where the weight function P(v) is a tempered distribution [50, 52]. Therefore, when con- 
sidering P-representable Gaussian density matrices, the weight function P(v) is simply a 
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FIG. 3: The •product of the width o~% with the power of the subfluctuant mode, for n — 10 1 (upper, 
blue), n — 10 2 (middle, red), and n = 10 3 (lower, black). For the horizontal axis, we used the 
non-linear scale x = arctan(ln(<5)) of Fig 2. For the vertical axis, we used the non-linear scale 
y = arctan(cr 2 ({Ilk, IT_k})). One clearly sees that for all n ^> 1, the product of Eq. (D16) stays 
between 1/4 and 1/2 as long as 8 < 1. One also verifies that it is asymptotically n-independcnt 
for all values of 5 -C n. Finally we recall that the power of the decaying mode coincides, to leading 
order, with ({II k , II_ k }), see Eq. (D17). 



Gaussian function. The generalization of this definition to two-mode (Gaussian) distribu- 
tions is straightforward: 



PP— repr 



d 2 v o9w 



7T 7T 



P(v,w) \v, k)(t>, k| (g> \w, —k)(w, — k| . 



(D20) 



The P-representation of the distribution characterized by Eqs. (15) can easily be cal- 
culated. From (D20), one calculates the expectation values Tr(nfc/5p- re pr) = n and 
Tr(dka_kAP-repr) = c, and one gets 



p ( v i w ) = ^7 ex P 



\ v \ i i * l \|2 

r7 \w — (cv n)\ 

n A' K ' n 



(D21) 



where A' 



n 



\c\ must be positive. The function P is well defined if, and only if, 



n > |c|, i.e. 5 > n/(n + 1). Hence, in the high occupation number limit, 

P - representability <=^ 5 > 1 . (D22) 
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This condition coincides with that for the non-existence of a sub-fluctuant variable, see Eq. 
(D14), as well as that for the vanishing of the off-diagonal matrix elements, see Eq. (D12). 
It also coincides with the condition obtained from the criterion of separability we shall 
discuss in the next subsection, see Eq. (D25). This proves that all quantum properties 
are lost if, and only if, S > 1. 

It is therefore of interest to further analyze this case and to determine how it is related 
to two-mode coherent states. Consider the limiting case A' — > + , i.e. 5 — > n/(n + 1) from 
above. In this limit P(v, w) of Eq. (D21) becomes 

8&) (w - V) . (D23) 

The presence of the two Dirac functions constraining the real and imaginary parts of w 
clearly expresses the critical character of this distribution. 

One can now notice that this distribution corresponds to the simplified distribution 
of Eq. (47). More precisely, the equality of the entropies, see Appendix D.7, suggest 
that (D23) is the quantum counterpart of the classical effective classical distribution Eqs. 
(26-28) obtained by setting the decaying mode to zero. The (badly named) minimally 
decohered distribution (44) is in fact slightly more decohered as can be understood from 
the replacement of the Dirac's 5^ (w — ^v*) by a Gaussian with a width equal to 1. This 
subtle difference is better appreciated when one realizes that the simplified distribution of 
Eq. (47) and that of (44), which are both diagonal in two-mode coherent states, belong 
to the same class of decohered matrices when considering the entropy. 

6. Separability of two-mode homogeneous distributions and violations of Bell's 

inequalities 

The non-separability is another important criterion for distinguishing distributions 
which cannot be viewed as classical stochastic distributions. Its physical relevance comes 
from the fact that the non-separability is a necessary property to have expectation values 
violating Bell's inequalities [9]. At the end of this subsection, we have included a brief 
discussion concerning the relationship between decoherence and the associated loss of the 
violation of Bell's inequalities. For more details we refer to [43]. 

A word of caution is in order. Unlike the P-representability the non-separability rests 
on a division of the 4 canonical variables into two subsets of 2 canonical variables. To define 
the two subsets, we use creation and destruction operators carrying a momentum k or — k. 
The criterion should therefore be called the (k, — k)-separability. The justification of the 
choice of the two subsets is presented after having explained the concept of separability. 

Following Werner [9], a distribution of two-modes k and — k is said to be separable if, 
and only if, it can be written as a convex sum of products of distributions, i.e. 

Pse P = Y.Pipl' '®/5-k. Pl>0, I> = 1- (D24) 
l i 

Notice that the definition (D24) is not the statement that a separable state is a product 
of its reduced density matrices (in which case one would obtain uncorrelated sub-systems: 
pi = for all I but one). 



Plim(v,w) = -exp 
n 



\v\ 
n 
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The key point is that all states of the form (D24) can be obtained by the following 
classical protocol: when a random generator produces the number / with probability pi, 
two distant observers performing separate measurements on the subsystems k and — k 
respectively, prepare them into the states p± k . Separable states are said classically cor- 
related because, by construction, their statistical properties can be interpreted classically. 
In particular they cannot violate Bell's inequalities. States which are not separable are 
called "entangled" (such as a two-mode squeezed vacuum states, or a singlet of 1/2 spins) 
and can only be produced by letting the two parts of the system interact. 

Although it may be very difficult to decipher the separability of a general two-mode 
state, there is a criterion for Gaussian states. We briefly outline the proof of this theorem, 
the complete proof can be found in [49]. First, one notices on (D24) that the separability 
is a property that is not affected by local unitary transformations XJ\ oc = <£> i/_k- These 
local transformations are of two types, rotations and squeezing which do not mix k and 
— k. Therefore, if p is a separable state, then all states p' = Ui oc pUi oc are also separable. 
Second, one notices that if p is a P-representable state (D20), it is separable. The non 
trivial part of the proof consists then to show that all Gaussian separable states can 
be constructed by applying local unitary transformations on a P-representable state. In 
brief, the criterion for separability for two-mode Gaussian states is a Gaussian state p is 
separable if, and only if there is a local unitary transformation Ui oc which brings it to a 
P-representable state. 

These local transformations are very constrained for homogeneous states because any 
local squeezing would break the homogeneity of the state (one would have either (d k ) / 
or (d^ k ) / 0, or both). Hence, Ui oc can only be a product of local rotations. These 
preserve the non-negativity of the Gaussian function P(v,w). As a result, separability 
and P-representability are equivalent properties of homogeneous Gaussian states. Given 
(D22), one concludes 

(k, — k) — separability of homogeneous states -4=>- S > 1 . (D25) 

Notice that the (k, — k)-separability and the ((pi, (^-separability are two independent 
properties of the state. The former implies the above inequality on the decoherence level. 
As shown in Appendix C, the latter is simply a consequence of homogeneity of the distri- 
bution. The independence of these two concepts can be understood by noticing that the 
separability defines an equivalent class of density matrices through local transformations 
only. In fact, the transformation from the operators with a definite momentum k to ((pi, 
02) is a global, 4 by 4, transformation, see (C5). 

We should now explain the selection of the two subsets. The choice of (dk, djj and 
(d_k> d^ k ), an d n °t (di, a\) and (d2, d!,), is dictated by the nature of the interactions 
amongst modes. When keeping non- linear terms, the Lagrangian density is a sum of 
products of fields evaluated at the same point. This leads to terms of the form (fikfiptp-k-p 
which are governed by the Fourier transform of the field (p(r],x). These non-linear terms 
give rise to momentum exchanges which cannot be described by di or d2 separately. (It 
should be noticed that these hermitian operators obey unusual commutation relations with 
the momentum operator P: [P, ai 2] = = Rkd2 > i). 

Notice that contrary to non-linear interactions, bilinear couplings do not permit to 
dismiss (pi, (f>2- Indeed, consider a bilinear interaction with a random field by adding 
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a coupling term Jd 3 x (j)(r), x) *&(r], x) where ^(rj, x) is a real stochastic field. Then the 
preferred basis are still the real and imaginary parts of 0(ry,x), see [21]. To verify this it 
suffices to decompose the interaction term in Fourier transform, use the reality condition of 
the stochastic field ^(77, x) to get fd 3 k 0i (77, k) (77, k) + ^(v, k) ^2(7?, k) , and then use 
the homogeneous and isotropic character of the moments of ^(r], x). In this we recover the 
fact that homogeneity implies that all Gaussian density matrices are (c/>i, <^>2)-factorizable. 

It is nevertheless possible to consider observables which cannot be written in terms 
of d\ or 02 separately without introducing explicitely non-linear interactions. To this 
end one should exploit the cosmic variance, namely the fact we do not observe the mean 
value, but rather some realization of the ensemble. (We refer to Chap. V of [33] for 
the procedure to extract some realization from the distribution, to Chap. VI for the 
calculation of observables associated with a given realization, and to the references of 
this work for the relationship between the role played by the interaction Hamiltonian in 
non-linear theories and the role of the projector introduced to extract some configurations 
from the mean.) In particular, upon specifying that a localized classical field configuration 
possessing a definite velocity (e.g. a localized wave packet made with running waves) has 
been realized/detected, the conditional expectation value of (/>(r/,x) describes a pair of 
wave packets moving away from each other and which are spatially well separated. These 
wave packets localized both in x and k space are built with and a_k separately. Hence 
they cannot be written in terms of a\ or a>2 separately. 

Loss of violations of Bell inequalities 
In [42], the authors point out that the projectors on a given field configuration can violate 
Bell's inequalities when evaluated in two-mode squeezed states. In [43], we extend their 
analysis by replacing the pure states by the partially decohered distributions of Eq. (D3) . 
We plot in Fig. 4 the combination (n + 1)C of Eq. (11) of [43] as a function of the 
square field amplitude for n = 100 and for three different values of 5. The points over the 
horizontal line (n + 1)C = 2 give combinations of operators violating Bell's inequalities. 
We observe that for 5 > 0.1, there is no more violations. However, the combination C 
does not provide the operator which maximizes the violation. Hence, there might still be 
operators which violate Bell's inequalities for 0.1 < 5 < 1. 



7. Entropy of partially decohered distributions 

To characterize the entropy of the general distribution (D3), we write 5 = An' 1 , A > 
0, 1 > 7 > 0. As we shall see, in the high occupation number limit, 7 is tightly constrained. 
The thermal occupation number defined in Eq. (51) and the entropy respectively scale 

as 

n = ^An^' 2 + 0(n-( 1+ T)/ 2 ) , (D26) 

and 

4?-k = (! + 7) 2r + ( 2 + lnA ) + 0{n-^' 2 ) . (D27) 

We have used the squeezing parameter n = sh 2 (r). From this analysis we see that the 
uncertainty in the definition of the quantum distributions which gives rise to the entropy 
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FIG. 4: The loss of violation as decoherence increases. The function is normalized to its maximal 
classical value = 2. The variable x gives the square held amplitude. The occupation number is 
n = 100, and the three values of S are (red, upper), the pure state, 0.05 (blue, middle), and 0.1 
(brown, lower), the border regime. It is also worth to point out that the function (n + l)C(x) is 
independent of n in the large n limit as what is found in Figs. 3, 5, and 6. 



of the effective distribution of sine functions of (28) is very limited. (The entropy of that 
distribution is Sk,-k = 2r, see Eq. (53) and the discussion in the following paragraphs.) 
Indeed, the linear dependence of the entropy in r implies that the distributions with 
entropy = 2r all have 

7 « 7^ • (D28) 
inn 

When refining the analysis by taking into account the term in In A, and requiring that 
15k, -k — 2(r + 1)| < B, one gets | lnS\ ~ 6 — 1 < B, independently of n. In the large n 
limit, this fixes the value of the coherence term c with a precision B/n, thereby tightly 
restricting the quantum density matrices which can be put in correspondence with Eq. 
(28). 

Nevertheless, the criterium Sk, -k = 2r + O(l) does not discriminate between the 
distribution of Eq. (44) (5 = 2) and that of Eq. (47) (5 = 1). Hence the distribution 
(44) should be conceived as minimal only as being a member of the equivalent class of 
minimally decohered matrices, i.e. those which have lost their quantum properties (5 > 1) 
and which possess an entropy equal to 2r. 



8. Partial randomization of the phase 



To complete the analysis we now assume that the phase of the squeezed vacuum state 
is randomly scattered. The density matrix pi n is then replaced by the statistical mixture 

f-2n 

p{r,<Po)= dipf(ip;ip ,a)p in (r,ip), (D29) 
J o 
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FIG. 5: Entropy as a function degree of coherence for three values of n. The function is normalized 
to it maximal value. The occupation number is n = 10 2 (lower, blue), n = 10 3 (middle, red), and 
n = 10 5 (upper, black). On the left panel, as a function of the ratio ( = |c|/(n + 1/2), one clearly 
see the sensitivity for £ — > 1 _ . On the right panel, as a function of x = arctan(ln(£)), one verifies 
that the critical value S ma x/2 corresponds to 5 = 1 for all n. 

where each pure state pi n (r, if) is parameterized by the squeezing parameters r and ip, 
and f(ip; ipo, a) is the PDF of ip centered at ipo with a spread a. We assume that the 
mean phase shift vanishes. Therefore the mean value of p> is the phase of the original 
distribution (12), namely (fo = tt + 2kr/ r . 

Since all the states in the ensemble (D29) have the same value of r, the expectation 
value of the occupation number is still (a^o-tk) = sh 2 r = n. However the effective 
coherence term is reduced since it is affected by averaging over the phase : 

{{c))y = Tr(a k a_ k/ S(r,v?o)) = y/n(n + l) / d(pf((p)e^. (D30) 

Jo 

Let us assume that the function f(<p) is sharply peaked. Then we can extend the 
bounds of integration to infinity. Let us choose for instance a Gaussian PDF, f{<p) = 
(2vrcr 2 )- 1 / 2 exp (-(ip - ip ) 2 /2a 2 ) . One gets 

{{c)) v = ^n(n + l)e^e- ff2 / 2 . (D31) 

Hence, we have that |((c}}^| = \c sq (r)\e~ a2 / 2 < \c sq (r)\, i.e. the norm of the effective 
coherence term is smaller, as expected. A value 5 cr = 1 corresponds to an extremely 
narrow distribution: a 2 r = l/(n + 1) = 4(A;77 r ) 4 ~ 10~ 100 . It should be noticed that these 
changes affecting the coherence are much smaller than those obtained when considering 
the trans-Planckian problem, see [27, 39, 40], which concerned modifications of the power 
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FIG. 6: Parametric curves We plot the product (j^({IIk, n_k}) versus the entropy. The latter 
(horizontal axis) has been normalized to its maximum value S ma x- The occupation number is 
n = 10 2 (upper, blue), n — 10 3 (middle, red), and n = 10 5 (lower, black). On the left panel, 
8 goes from to n and the product cr 2 ({llk, n_k}) has been normalized to its maximum value 
n. On the right panel, we show a zoom on the region S < 1, and the product is not normalized. 
One sees the two regimes. First as long as the entropy has not reached half its maximum value, 
the product stays nearly constant. On the contrary, for S > S max /2, the product explodes and 
reaches n. One can also admire the fact that the curves meet near (1/2, 1/2) irrespectively of the 
value of n. Moreover in the limit n — > oo the curve becomes an Heavyside function, thereby cleary 
separating the two regimes: the first one where quantum properties are progressively lost, and the 
second one where the power of the decaying mode linearly grows with 6. 



n(k). 
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